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ABSTRACT 

Detection of caustic crossings of binary-lens gravitational microlensing events is im- 
portant because by detecting them one can obtain useful information both about the 
lens and source star. In this paper, we compute the distribution of the intervals be- 
tween two successive caustic crossings, f(t cc ), for Galactic bulge binary-lens events to 
investigate the observational strategy for the optimal detection and resolution of caus- 
tic crossings. From this computation, we find that the distribution is highly skewed 
toward short t cc and peaks at t cc ~ 1.5 days. For the maximal detection of caustic 
crossings, therefore, prompt initiation of followup observations for intensive monitor- 
ing of events will be important. We estimate that under the strategy of the current 
followup observations with a second caustic-crossing preparation time of ~ 2 days, 
the fraction of events with resolvable caustic crossing is ~ 80%. We find that if the 
followup observations can be initiated within 1 day after the first caustic crossing by 
adopting more aggressive observational strategies, the detection rate can be improved 
into ~ 90%. 
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INTRODUCTION 



Several groups are searching for massive astronomical com- 
pact objects (MACHOs) by monitoring light variations 
of stars caused by gravitational microlensing towards the 
Galactic bulge and the Magellanic Cloud fields (EROS: 
Aubourg et al. 1993; MACHO: Alcock et al. 1993; OGLE: 
Udalski et al. 1993; DUO: Alard & Guibert 1997). With their 
efforts, more than 400 candidate events have been detected 
to date. 

Among these events, a considerable number of events 
are believed to be caused by binary lenses (Udalski et al. 
1994a; Dominik & Hirshfeld 1994, 1996; Mao & Di Stefano 
1995; Alard, Mao, & Guibert 1995; Bennett et al. 1996; Rhie 
& Bennett 1996; Alcock et al. 1999b; Afonso et al. 2000). 
The light curve of a binary-lens event differs from that of a 
single- lens event. The most distinctive feature of a binary- 
lens event light curve occurs when a source star crosses a 
lens caustic. When the source star is located on the caustic, 
the amplification becomes very large. Therefore, the light 
curve of a caustic-crossing binary-lens event is characterized 
by a sharp spike. On the other hand, if a binary-lens event 
does not involve a caustic crossing, the resulting light curve 
exhibits a relatively small deviation from that of a single- 
lens event (Mao & Paczyriski 1991; Dominik 1998). 



Detection of caustic-crossing binary-lens events is im- 
portant because they can provide useful information both 
about the lens and source star. First, a caustic-crossing event 
provides an opportunity to measure how long it takes for the 
caustic line to transit the face of the source star. By using 
the source crossing time along with an independent deter- 
mination of the source star size, one can determine the lens 
proper motion with respect to the source star (Afonso et al. 
1998; Albrow et al. 1999a; Alcock et al. 1999a). With the 
determined lens proper motion, one can significantly better 
constrain the physical parameters of the lens (Gould 1994; 
Nemiroff & Wickramasinghe 1994; Witt & Mao 1994; Peng 
1997). Second, a caustic-crossing event can also be used to 
determine the surface structure of the source star such as 
the radial brightness profiles (Gaudi & Gould 1999; Albrow 
et al. 1999b, 2000; Afonso et al. 2000) and spots (Han et 
al. 1999). To obtain these useful information from the light 
curve of a binary-lens event, the event should be monitored 
with high time resolution. However, under the nightly mon- 
itoring observational strategy of the current lensing experi- 
ments, it is difficult to construct light curves with resolution 
high enough to resolve caustics. 

However, the caustic crossing can be resolved with the 
help of alert system based on real time observations (Alcock 
et al. 1996, 1997; Udalski et al. 1994b) and subsequent high 
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time resolution followup observations (GMAN: Pratt et al. 
1996; PLANET: Albrow et al. 1998; MPS: Rhie et al. 1998). 
Since the caustics of a binary-lens event form a closed curve, 
the source star of a caustic-crossing event crosses the caus- 
tic line at least twice. Although the first caustic crossing is 
unlikely to be resolved due to its short duration, it can be 
inferred from the enhanced amplification. Then, if followup 
observations can be prepared before the second caustic cross- 
ing, dense enough sampling through the second caustic will 
be possible. Therefore, it is important to estimate the dis- 
tribution of the intervals between caustic crossings (caustic- 
crossing intervals, t cc ) for the optimal detection and resolu- 
tion of caustic crossings. 

The distribution of caustic-crossing intervals, f(t cc ), 
was computed by Honma (1999) for events expected towards 
the Magellanic Cloud field. His intention for the computa- 
tion of f(t C c) was to demonstrate the detection bias against 
events with short i cc , for which the second caustic cross- 
ings are more likely to be missed. For this purpose, it was 
enough to determine f(t cc ) based on the typical values of the 
lens mass and the Einstein ring radius crossing time (Ein- 
stein time scale, t^)- However, for the investigation of the 
optimal observational strategy to detect and resolve caus- 
tic crossings, it is required to determine /(£ C c) based on the 
distribution of Einstein time scales expected from detailed 
models of the physical and dynamical distributions of lens 
matter and binary mass function. In addition, since majority 
of caustic-crossing binary-lens events have been detected to- 
wards the Galactic bulge (Alcock 1999b), constructing f(t cc ) 
for these events is desirable. 

In this paper, we compute the distribution of the in- 
tervals between two successive caustic crossings, f(t cc ), for 
Galactic bulge binary-lens events to investigate the observa- 
tional strategy for the optimal detection and resolution of 
caustic crossings. From this computation, we find that the 
distribution is highly skewed toward short t cc and peaks at 
t cc ~ 1.5 days. For the maximal detection of caustic cross- 
ings, therefore, prompt initiation of followup observations 
for intensive monitoring of events will be important. We es- 
timate that under the strategy of the current followup ob- 
servations with a second caustic-crossing preparation time 
of ~ 2 days, the fraction of events with resolvable caustic 
crossing is ~ 80%. We find that if the followup observations 
can be initiated within 1 day after the first caustic crossing 
by adopting more aggressive observational strategies, the 
detection rate can be improved into ~ 90%. 



2 CAUSTIC-CROSSING INTERVALS 
2.1 Binary-lens Events 

When the lengths are normalized to the combined Einstein 
ring radius, the lens equation in complex notations for a 
binary-lens system is represented by 

TTll 7T12 



c = 



Zl 



+ 



«2 



(1) 



where mi and m,2 are the mass fractions of individual lenses 
(and thus mi + mi = 1), z\ and Z2 are the positions of the 
lenses, £ = £ + ir\ and z = x + iy are the positions of the 
source and images, and z denotes the complex conjugate of 
z (Witt 1990). The combined Einstein ring radius is related 



0.5 - 




Figure 1. Caustic crossings for various source star trajectories 
(thin solid straight lines) for an example event caused by a binary 
system with a separation I = 1.0 and a mass ratio q = 1.0. The 
caustics are represented by a thick sold line. The location of the 
lenses (mi and 1112) are marked by dots and the circle drawn by 
a dashed line represents the combined Einstein ring. 



to the total mass M of the binary system and its location 
between the observer and the source star by 



= / AGM D ol D lt 
~ V c 2 D os 



1/2 



(2) 



where D a i, Di 3 , and D a3 are the separations between the ob- 
server, lens, and source star. Since the individual microlens- 
ing images cannot be resolved due to the small separations 
between them, the amplification of the binary-lens event is 
given by the sum of the amplifications of the individual im- 
ages, Ai, which are given by the Jacobian of the transfor- 
mation (1) evaluated at the image position, i.e. 



At = 



( 1 



V |det J\ 



det J 



dz dz 



£3- (3) 



The caustic refers to the source position on which the am- 
plification becomes infinity, i.e. det J = 0.Q The set of caus- 
tics form a closed curve, called cautsics. The caustics take 
various shape and size depending on the binary mass ratio 
q — mi/m,2 and the projected binary separation I normal- 
ized by the angular Einstein ring radius #e = tb/D i. In 
Figure 1, we present the caustics (a thick solid curve) of 
an example binary- lens event with I = 1.0 and q = 1.0. 
In the figure, the locations of the lenses (mi and 7712) are 
marked by dots and the circle with its center at the cen- 
ter of mass of the binary system drawn by a dashed line 
represents the combined Einstein ring. For more example 



* The actually observed amplification is finite because the source 
star is not a point source. The amplification for an extended 
source is the weighted mean of the amplification factors over the 
surface of the source, leading to a finite value (Schneider & Weiss 
1986). 
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caustics of binary-lens events with various values of t and q, 
see Figure 1 of Han (1999). 



2.2 Distribution of Caustic-Crossing Separations 

The caustic-crossing interval is proportional to the separa- 
tion between two successive caustic-crossing points (caustic- 
crossing separation, 9 CC ). Due to the dependency of the 
shape and size of the caustics on the binary separation and 
the mass ratio, the caustic-crossing separation and the re- 
sulting caustic-crossing interval depend on the values of £ 
and q. In addition, since t cc is scaled by the Einstein time 
scale, which depends on the physical lens parameters of the 
mass, the location, and the lens-source transverse speed v 
by 



q = 1.0 q=0.B q=0.6 q 0.4 q 0.2 



t H = / 4GM DpiDis 

V \ C 2 V 2 Dos 



1/2 



(4) 



the caustic-crossing interval depends also on these parame- 
ters. Therefore, for the determination of /(tcc), one should 
consider binary-lens events which are expected for all com- 
binations of £ and q over the entire range of the physical lens 
parameters. 

To determine f(t cc ), we first determine the distribu- 
tion of the caustic-crossing separations f(9 cc ;£,q), which is 
expected for events caused by a binary lens with £ and q. 
By defining a binary-lens event as a close lens-source en- 
counter within the combined Einstein ringjj], we produce a 
large number of source trajectories that pass through the 
combined Einstein ring. The source star trajectory orienta- 
tions with respect to the projected binary axis, ip, and the 
impact parameter, /3, are randomly selected in the ranges 
of < ip < 2-7T and < j3 < 1, respectively. If the source 
star trajectory crosses the lens caustics, we then find the 
crossing points and measure the separation 9 CC in units of 
#e- In most cases, the source crosses the caustics twice, e.g. 
the three trajectories numbered by 1, 2, and 3 in Figure 1. 
But due to the concavity of the caustics, crossings can occur 
more than twice, e.g. the trajectory numbered by 4 in Fig- 
ure 1. For these cases, we measure 9 CC between every pair of 
two successive crossing points. 

In Figure 2, we present the distributions f(9 cc ;£,q) for 
events caused by binaries with various projected separations 
and mass ratios. The individual distributions are arbitrar- 
ily normalized, but they are relatively scaled so that the 
area under each distribution is proportional to the caustic- 
crossing probability. From the figure, one finds that the dis- 
tribution f(9 cc ;£, q) depends strongly, both in the scale (i.e. 
the caustic-crossing probability) and the shape, on the bi- 
nary separation. In the scale, the caustic probability be- 
comes maximum for binary- lens events with £ ~ 1.5. This 
is because caustics at around this binary separation form 
the largest curve. With the same reason, the distribution 
peaks at relatively large caustic-crossing separations, e.g. at 



t Some binaries form their caustics outside their combined Ein- 
stein ring. However, the size of the outer caustics is usually very 
small compared with the combined Einstein ring, implying that 
the probability of the outer caustic crossing will be very small. 
Therefore, determination of f(9 cc ;£,q) will not be seriously af- 
fected by our definition of a binary-lens event. 




Figure 2. The distributions of the normalized caustic-crossing 
separations for events caused by binary systems with various sep- 
arations £ (also normalized by 0e) and mass ratios q. The indi- 
vidual distributions are arbitrarily normalized, but they are rel- 
atively scaled so that the area under each distribution is propor- 
tional to the caustic-crossing probability. 



around 9 CC = 0.4 for binary-lens events with £ ~ 1.5. As the 
separation becomes smaller or larger, the caustics becomes 
smaller, making the caustic-crossing probability smaller and 
the distribution peaks at small value of 9 CC . On the other 
hand, the dependency of f(9 cc ;£,q) on the binary mass ra- 
tio both in the scale and the shape is relatively weak (Gaudi 
& Gould 1999). 

With the obtained distribution of f(9 cc ;£,q), we then 
obtain the combined distribution of caustic-crossing separa- 
tions by 



f(0cc) 



dqf(q) / dbf(9 cc ;£,q)f(b)S(9 cc 



(5) 



where f(q) and f(b) respectively are the distributions of the 
binary mass ratios and the intrinsic binary separations b in 
units of 6*e, and S denote the delta function. The projected 
binary separation £ is obtained from the intrinsic separation 
b by assuming that the intrinsic binary separation vector b is 
randomly oriented onto the sky. The notation 9' cc represents 
the caustic-crossing separation obtained under the assump- 
tion of random orientation of b. For the distribution f(q), 
we adopt the distribution determined by Trimble (1990). 



2.3 Distribution of Caustic-Crossing Intervals 

Once the distribution of caustic-crossing separations f(9 cc ) 
is obtained, the distribution of the caustic-crossing intervals 
expected for events caused by binaries with a constant mass 
M is obtained by 



/(4) = 



dD os p(D os ) 



( D iDi s \ 1/2 



© 0000 RAS, MNRAS 000, 



4 Han, C, Park, S.-H., & Lee, Y.-S. 



dv z vf(v y ,v z ) / d9 C cf{0 C c)8(t' cc -6 cc t E ) 



dv,, 



(6) 

where p(D os ) and p(D a i) are the mass density of source stars 
and lenses along the line of sight toward the Galactic bulge, 
rfmax is the maximum extent of the source star distribution, 
(v y ,v z ) are the two components, parallel and normal to the 
Galactic plane, of the lens-source transverse velocity v, and 
f(v y ,v z ) represents their distribution. In the equation, the 
factors (DoiDis/Dos) 1 ^ 2 and v are included to weight f(t' cc ) 
by the cross section of the lens-source encounter, i.e. the Ein- 
stein ring radius r-E oc (D iDi a / Dos) 1 ^ 2 , and the transverse 
speed. 

For the matter distributions of the Galactic bulge, we 
adopt a 'revised COBE' model, which is based on a tri- 
axial COBE model (Dwek et al. 1995) except the central 
part of the bulge. For the inner ~ 600 pc of the bulges, we 
adopt a high central density Kent (1992) model which bet- 
ter matches with observations in this region. For the matter 
distributions of the Galactic disk, we adopt a double expo- 
nential disk model with vertical and radial scale heights of 
ha — 325 pc and h z = 3.5 kpc (Bahcall 1986). The trans- 
verse velocity distributions for both Galactic bulge and disk 
lenses are modelled by a Gaussian of the form 



f(vi) oc exp 



(Vi 



2(7, 



(7) 



where the means and the standard deviations for individ- 
ual transverse velocity components are (v y , a y ) — (220,30) 
km s _1 and (v z ,a z ) — (0,20) km s _1 for the disk lenses, and 
(i; y ,cr H ) = (220,93) km s" 1 and (v z , a z ) = (0,79) km s" 1 for 
the bulge lenses, respectively. For the detailed discussion 
both about the matter and the transverse velocity distribu- 
tions, see Han & Gould (1996). 

The expected distributions of f(t' cc ) for events caused 
by various values of the binary mass are presented in the 
upper panel of Figure 3. In the lower panel, we also present 
the integrated probability for events with t > t' cc , i.e. 
P(t > 4) = 1 - f* /(4)</ /(fcc)<*C- In each panel, 
to better show the distribution in the short caustic-crossing 
interval region, we expand the region and present in a sepa- 
rate small box. From the figure, one finds that as the binary- 
lens mass decreases, the distribution becomes narrower and 
peaks at shorter t cc . 

Binary-lens events are caused by Galactic binaries with 
various masses. Then the final distribution of the caustic- 
crossing intervals is obtained by convolving f(t' cc ) with the 
mass function of binary lenses <E>bi(M), i.e. 

/(*«)= / dMM 1/2 $ hi (M) / dt' cc f(t' cc )5(t cc ~M 1/2 t' cc ). 



(8) 

The factor M 1 ^ 2 in the first integrand is included to weight 
/(ice) by the cross-section of lens-source encounter (te oc 
M 1/2 ), while the same factor in the second integrand is 
included because the caustic crossing-crossing interval is 
scaled by the Einstein time scale (te oc M 1 ^ 2 ). 

The mass function of binary lenses is very uncertain 
due to the observational difficulties in identifying binaries 
by resolving their components combined with our ignorance 
about the mechanism of binary formation. Furthermore, 
since binary-lens microlensing events can be caused by bina- 
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Figure 3. Upper panel: The distributions of caustic-crossing 
intervals for Galactic bulge events caused by binaries with vari- 
ous constant masses. Lower panel: The integrated probability for 
caustic-crossing events with t > t' cc . In each panel, to better show 
the short t' cc region, we expand the region and presented in a 
separate small box. 



ries which are composed of dark components (either one or 
both), construction of $bi(A/) becomes even more compli- 
cated. Though many possibilities have been proposed, the 
scenarios for the formation of binaries can be classified into 
two categories. The first one is that binaries are formed dur- 
ing the star formation stage through fragmentation of a col- 
lapsing object (Norman & Wilson 1970; Boss 1988). In the 
other scenario, by contrast, both binary components form 
independently and then later combine through some sort 
of capture mechanism (Fabian, Pringle, & Rees 1975; Press 
& Teukolsky 1997; Ray, Kembhavi, & Anita 1987). Due to 
the difference in the binary formation mechanism, the mass 
functions expected from the two scanarios would be differ- 
ent each other. Under the formation environments of the 
former scenario, the most probable binary-lens mass func- 
tion would be similar to that of single lenses. On ther other 
hand, under the circumstances of the latter binary formation 
scenario, the two lens masses would be drawn independently 
from the same mass function, which is similar to the single- 
lens mass function. For the construction of /(tcc), therefore, 
we test both binary mass function models expected under 
the two cases of binary formation scenario. The first binary 
mass function is modeled by a power law with a mass cutoff, 
i.e. 



$bi(M) oc M" p e(m 



(9) 
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Figure 4. The model binary-lens mass functions that are used for 
the determination of the caustic-crossing interval distributions. 



where O is a heavy side step function. For the values of 
the power and the mass cutoff, we adopt p = 2.1 and 
m cut — 0.04 Mq, following the best-fitting values deter- 
mined from the distribution of Einstein time scales of 51 
Galactic bulge single- lens events by Han & Gould (1996). 
The second model of binary mass function is constructed by 
selecting masses of individual binary components from the 
same mass function in equation (9) and then combining the 
two masses. Since two masses are combined to yield a sin- 
gle binary mass, binaries following the latter mass function 
tend to be heavier than the binary lens masses following 
the former mass function. We distinguish the two models 
by calling them 'heavy' and 'light' binary mass functions. 
Note that since the mass cutoff of both the heavy and light 
binary mass functions is less than the hydrogen-burning 
limit of 0.08 Mq, dark component of lenses are included 
in these mass functions. In the third model, we test binary 
mass function which is constructed under the assumption 
that lenses are composed of only stars: 'stellar' mass func- 
tion. We construct the stellar mass function by adopting the 
mass function determined by Zoccali et al. (2000) from ob- 
servations of Galactic bulge stars by using the Hubble Space 
Telescope (HST) plus the Near Infrared Camera and Multi 
Object Spectrometer (NICOMOS). The stellar binary-lens 
mass function extends down to 0.15 Mq. In Figure 4, we 
present the model binary-lens mass functions. 

In Figure 5, we present the finally determined distribu- 
tion of caustic-crossing intervals and integrated probability 
with t > t cc . From the figure, one finds that the distribu- 
tions /(tec) are highly skewed toward short t cc regardless of 
the assumed mass functions. The peak of the distribution 
occurs at ~ 1.5 Mq although ther exist slight variations de- 
pending on the assumed binary-lens mass functions. As a 
result, majority (70%-90%) of binary-lens events expected 




20 40 
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Figure 5. The resulting distributions of caustic-crossing in- 
tervals (lower left panel) and the integrated probabilities with 
t > t cc for Galactic bulge binary-lens events determined for the 
two model binary-lens mass functions in the upper panel. 

to be detected towards the Galactic bulge will have caustic- 
crossing intervals shorter than t cc = 20 days. 



3 DISCUSSION 

In the current strategy, intensive followup observations of 
the second caustic crossing of a binary-lens event can be ini- 
tiated within ~ 2 days after the first crossing. The primary 
search teams alert on the first caustic crossing within 1 day, 
and the followup teams begin monitoring of the event im- 
mediately. For example, the PLANET group generally pro- 
duces template images the same day and is able to make 
rudimentary measurements of the event within a day af- 
ter the receipt of the alert from the primary search teams. 
Once these rudimentary measurements of the event are in 
place, one can determine a second crossing time without dif- 
ficulties. Therefore, the typical response time of the current 
followup observation teams is 2 days if the weather is good 
(A. Gould 1999, private communication). According to the 
determined distribution f(t cc ), then, the fraction of binary- 
lens events whose second caustic crossings can be monitored 
by the followup observations is ~ 80% although there exist 
slight difference depending on the assumed mass functions. 
Therefore, we find that for Galactic bulge events the detec- 
tion bias toward long time-scale events argued by Honma 
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(1999) will not be very strong. However, we note that his 
argument about the detection bias against short time-scale 
events is based on the assumption that the current prepara- 
tion time for intensive followup observations is greater than 
7 days. We are not certain why his adopted preparation time 
is that long. However, if the same preparation time of 2 days 
is adopted, the fraction of LMC binary-lens events that can 
be intensively monitored by followup observations would be 
similar to that of Galactic bulge events. 

However, the expected most frequent caustic-crossing 
interval of ~ 1.5 day is shorter than the current preparation 
time of the followup observations. Then, if the preparation 
time is shortened by adopting more aggressive observational 
strategy, one can detect and resolve caustic crossings for 
more binary-lens events. With the help of automatized pro- 
cess, the MACHO group (Alcock et al. 1997) can finish data 
reduction within 6 hrs of image acquisition. With vigorous 
efforts to minimize the preparation time, therefore, initiation 
of followup observations within 1 day after the first caustic 
crossing will be possible. We find that if the followup obser- 
vations can be initiated within 1 day after the first caustic 
crossing, the detection rate can be improved into ~ 90%. 
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